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METHOD OF ANALYSIS FOR COMPRESSIBLE FLOW THROUGH MIXED-FLOW CENTRIFUGAL 
IMPELLERS OF ARBITRARY DESIGN ! 


By Josera T. HAMRICK, AMBROSE GINSBURG, and WALTER M. OSBORN 


SUMMARY 


A method is presented for analysis of the compressible flow 
between the hub and the shroud of mixed-flow impellers of 
arbitrary design. Axial symmetry was assumed, but the forces 
in the meridional (hub to shroud) plane, which are derived from 
tangential pressure gradients, were taken into account. 

The method was applied to an experimental mixed-flow im- 
peller. The analysis of the flow in the meridional plane of 
the impeller showed that the rotational forces, the blade curra- 
ture, and the hub-shroud profile can introduce severe velocity 
yradients along the hub and the shroud surfaces. Choked 
fow at the impeller inlet as determined by the analysis was 
verified by experimental results. 


INTRODUCTION 


Because of the complexity of the numerical solution for a 
three-dimensional analysis of flow through a centrifugal com- 
pressor, most treatments of the problem have been concerned 
with a two-dimensional analysis of the flow between the 
pressure side of one impeller blade and the suction side of 
the preceding blade. Relaxation methods in which it was 
assumed that no velocity variation existed between the hub 
and the shroud have been used (reference 1) to give an 
analysis of the flow between blades. In the conventional 
mixed-flow centrifugal impeller, however, the nonsymmetri- 
eal bub-shroud profile, the nonuniform energy addition 
across the passage, and the blade-angle variation due to close 
adherence to radial blade elements may introduce velocity 
gradients between the hub and the shroud comparable with 
those between the blades. For mixed-flow compressors, 
therefore. an analysis of the flow in the meridional plane is 
required. (The meridional plane is one that passes through 
and is parallel to the axis of rotation.) 

In order to investigate the isentropic compressible flow in 
the meridional plane, a fow equation based on equilibrium, 
continuity, and. general energy considerations was derived 
at the NACA Lewis laboratory during 1949-50 and is 
presented herein, Axial symmetry was assumed, but the 
forces that result from tangential pressure gradients were 
taken into account. The numerical solution involves the 
adaptation of stream-filament theory and the numerical 
integration of all equations. 

The applicability of the method for predicting the flow 
within the impeller was investigated by analyzing the flow 


through a mixed-flow centrifugal impeller for which experi- 
mental data are available (reference 2). Flow streamlines, 
pressure distributions, and velocity distributions relative to 
the impeller for a given set of operating conditions are 
presented. Theoretical shroud pressures and theoretical 
maximum weight flows over a range of speeds are compared 
with those experimentally obtained. 


ANALYSIS 


An analysis is made in the meridional plane of the mixed- 
flow impeller. Axial symmetry is assumed, but forces in the 
meridional plane that are derived from tangential pressure 
gradients are taken into account. By the use of the methods 
presented in this analysis, the isentropic compressible flow 
in the meridional plane may be obtained and shown in 
terms of the flow streamlines and the velocity and pressure 
distributions. 






-~-Normal line- 
"= Fluid particle. - 


“~- Streamline --~---- 








rotation 






Conic surface 
generated by 
tangent tine-4 





(5) 


(a) Impeller profile in meridional plane. 
(b) Intersection of developed cone surface. 


FIGURE I.— View of impeller showing orientation of blade angle and forces acting on rotating 
fuid particle in meridiona] plane. 


t Supersedes NACA FN 2165, “Method of Analysis for Compressible Fiow Through Mixed-Flow Centrifugal Impellers of Arbitrary Design” by Joseph T. Hamrick, Ambrose Ginsburg, 


and Walter M. Osborn, 1950. 
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Geometrical quantities —A complete geometrical. descrip- 
tion of the blade is necessary for a numerical solution. 
Geometrical quantities are defined in appendix A and shown 
in figures 1 and 2. The blade angle £ is oriented in the 
following manner: In figure l(a), a line OB, which intersects 
the axis of rotation at O, is drawn tangent to the pomt P on 
a streamline in the meridional plane. If the tangent line 
OR is rotated about the center line of impeller rotation 0-0’, 
a conic surface of revolution that intersects the impeller 
blades is formed. This intersection is shown on a developed 
surface of the cone in figure 1(b). At the point of tangency 
to the streamline (point P), the blade mean line forms the 
angle —8 with the tangent line. 


Limitations and assumptions.—-The fluid is assumed to be 
nonviscous but compressible. All compression is by the 
adiabatic, isentropic, steady-flow process. The velocity is 
assumed uniform from blade to blade in the direction of 
rotation (axial symmetry) and the average angle of flow 
relative to the impeller is assumed the same as the blade 
angle 8. This assumption is considered valid where channel 
flow prevails, but is questionable at the inlet and the outlet. 
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(ec) Orientation of e with respect to 8. (Dashed lines Jie in meridional plane.) 


FIGURE 2.— Quantities used in development of force equations. 


Method of analysis.—-The method 1s based on one proposed 
by Flügel for analysis of flow in the blade-to-blade direction 
and discussed by Stodola in reference 3. The forces acting 
on a particle along the streamline path are considered. The 
component of the blade surface force normal to the meri- 
dional streamline, the centrifugal forces resulting from the 


curvilinear path of the particle, and the force due to rotation 

about the axis of the impeller are equated to the change in 

pressure across the particle in the hub-to-shroud direction. 

The resulting force equation is combined with one based on 
the general energy equation to give an equation involving 

the relative velocity g and a number of variables dependent 

upon the geometry of the impeller. The variables that lend. 
themselves to graphical treatment are lumped into two para- 

meters giving a first-order differential equation linear in g. 

This equation is used in conjunction with the continuity 

equation to obtain the velocity distribution. Streamlines 

between the hub and the shroud similar to those shown in 

figure 1 (a) are drawn in by estimate. When the differential 

equation is solved for q and a velocity is assumed at the hub, 

the velocity distribution from hub to shroud may be obtained 

by numerical integration. By means of the continuity equa- 
tion, the weight flow corresponding to the assumed hub 

velocity is obtained. Several solutions may be necessary for 
obtaining the desired weight flow. The streamlines may be 
rearranged to match the new values of velocity and density 
and the complete analysis repeated. 


Force equation.—For the forces acting on the fluid 
particle in the n direction on the meridional plane (fig. 1 (a)), 
the following equation is obtained (appendix B): 


ldp_ 
pdn = r 


? (artu)? Aa F, (1) 


From figure 1 it can be seen that 


v=g cos Bp 
and 
u=g sin f 


In terms of g, equation (1) becomes 


ldp (— cos? 8 4 Sin’ sin? 8 dr 
pan 4 r. r an 


d A a 
inte r+ oqsin pF (2) 


Energy addition——From general energy considerations, 
the isentropic adiabatic energy addition is 


; [wry (wry t ud wr) = fa z (ar Hy +o: — H; 


where the prerotation term \=7;(wr;-+t,) 


aE 2 emma (3) 





From figure 1, 
gmurt py? 


and with the assumption that H, is constant for all stream- 
lines, 


wr; 
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For an isentropic process 


T Pee ee, 
-r= | > Pp 
Therefore 
aL 
ae dp= z; luri rA a gA aAA M) (4) 


When differentiated with respect to n, equation (4) becomes 


ldp_, dr dg dù 
aa an oda an (5) 


General velocity equation.—Equations (2) and (5) are 
combined to give the following equation: 


dq__ (cos? ms dr F, od 
dn ( Pe en sing Fe q gan (6) 


For straight blades lying on the meridional plane with no 
prerotation of the fluid at the inlet, the terms fa and 2 e 
drop out; a first-order differential equation linear in q oe 
remains. In order to arrive at a solution for curved-blade 
impellers and for prerotation at the inlet, it may be desirable 
to solve for g and to neglect these terms in the first solution. 
The values of ¢ thus found may be used in the force and pre- 
rotation terms for the second solution and may be treated in 
the same manner as those terms dependent upon the geometry 
of the impeller. 


Velocity equation for impeller of arbitrary blade shape 
having radial blade elements with no prerotation.—The 
blade-force component F;, normal to the meridional stream- 
lines is evaluated in appendix C as 


a a eee 
ri f tan æ (20 sin ata + sin’ sin @ tan æ 


The velocity equation for a radial-bladed impeller without 
prerotation becomes 


dg___fcos*p sin*ß dr _ sin?’ )- 
dn = (15E re å r da r ee 


Qu sin 8 or sin Sana (20 gi at) (7) 


Variables dependent upon the impeller geometry or upon 
values obtained from a previous solution are represented by 
parameters a and b where 


_ cos? sin? 8 dr sin?8 . 
cos’ A __ sin’ g B sin a tan æ 
re r dn r 





and 
b=2w sin 8 ot sin 8 tan æ (2 sin atge) 


Equation (7) becomes 


ih PETR 
aad b (8) 


When equation (8) is solved for g, 


yao a Ce re an] (9) 


In solving for q, the term du/dm in parameter b is neglected 
for the initial solution in order to obtain estimated values of 
du/dm for the succeeding solution. 

Blade surface velocity distribution.—If{ a linear pressure 
gradient from blade to blade is assumed, the blade surface 
velocity distribution may be approximated for a finite number 
of blades by the method given in appendix D. 

Continuity equation.—The equation for continuity of 
flow through the hub-to-shroud passage may be written 


w=2r | rfoga cos 8 dn (10) 


where f is a factor that takes into consideration the area 
taken up by the blade cross section. The value of the 
density p is determined from the energy equation (equation 
(3}) and the relations 
i)" 
on i H, 


paris OROS A 


NUMERICAL EXAMPLE 
IMPELLER EXAMPLE 





(11) 


A numerical example has been computed for an 18-blade 
mixed-flow centrifugal impeller, the performance of which is 
described in reference 2. The impeller was designed for a 
constant cross-sectional flow area from inlet to outlet. The 
blade curvature is such that the intersection of the blade 
with any circular cylinder having the same center line as the 
impeller produces a parabola on the developed cylinder sur- 
face. Because of adherence to radial blade elements, the 
blade angle 8 increases from hub to shroud. 

The operating conditions for the impeller were as follows: 


Equivalent impeller tip speed, U/4/@, ft/see..._...___-_-_____ 133i 
Ratio of specific heats, y__------.--.-----_--~--.----+--~---- 1.4 
Prerotation At Mit- c 225220 oie eee oe ke oes 0 
Impeller equivalent weight flow, 14/8, Ib/sec_..._...._.. 8.743 


PROCEDURES FOR NUMERICAL SOLUTION 


The meridional plane of the impeller together with the 
estimated streamlines and resulting normals is shown in figure 
l(a). Estimated streamlines are drawn in between the hub 
and the shroud; equal mass flow is assumed between the 
annular stream tubes. Lines normal to the streamlines are 
drawn in along the flow path. 

Radius of curvature.—The radius of curvature of the 
streamline at any point in the meridional plane r, is obtained 


from. 
5 _ G+ 
— r” 


The slope of each streamline r’ at each normal v is plotted 
against axial distance in order to obtain the change in slope r”. 
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Numerical integration—The velocity q@ is found by 
numerically integrating equation (9) 


neigh “(or fre! an] (9) 


The velocity g is the estimated velocity at the hub. 
When a is plotted against n (fig. 3(a)) where n is the dis- 
tance along the normal from hub to shroud, the integral 


ii a dn can be evaluated by the trapezoidal rule. <A plot of 


aoe against n can then be obtained (fig. 3(b)) and 
ey ied in & similar manner. As previously stated, the 
term du/dm in parameter 6 is neglected in the first solution. 
The evaluations of these integrals allow the determmation of 
the velocity g at any point n. 

The weight flow at any normal may be obtained from 
equation (10) 


w=2r ["rfeea cos 8 dn (10) 


When a plot of rfpgg cos 8 against n is used (fig. 3(c)), the 


integral f : rfpgg cos 8 dn can be evaluated by the trapezoidal 


rule. If the weight flow does not equal the designated 
weight flow for the passage, a new value of the velocity go 
is assumed and the solution is repeated. 

In order to satisfy continuity for each stream tube, new 
streamlines are drawn based on the velocities of the preceding 
solution. Velocity q (fig. 4(a)) and density pg (fig. 4(b)) 
are plotted against n. The radii at which the new stream- 
lines cross the normals are obtained as follows: From 


Parameter a, ft 
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A fry fa 
COS @ 
and 
Aw=2arfpgq cos f An 


AW COS a 
— 2 
hg Va FF peq cos B 


where cos a, f, p, g, and cos 8 are average values between 
streamlines based on estimated values of fa} After the 
new streamlines are determined, the normals are drawn. 
The graphical solution is then repeated for succecding scts 
of new streamlines until continuity is satisfied. 


ACCURACY 


The accuracy of a solution depends to a large extent upon 
the accuracy with which the streamline curvature is deter- 
mined. Measurement of other geometrical quantities, of 
course, enters into the accuracy, but their measurement is 
more precise. 


The sudden decrease in flow area (as much as 25 percent 
at the hub) from & point just ahead of the blades at the inlet 
to a point just inside the blades distorts the streamline pic- 
ture. A tedious and detailed solution is necessary for a true 
picture of the distortion in this region. Also, the assumption 
that the flow angle is the same as that of the blades at the 
inlet is questionable. A lengthy analysis at the inlet was 
therefore considered unjustifiable. 

The cross-sectional area for the clearance is approximately 
4 percent near the inlet and affects the maximum flow. 
In the solution, the clearance area was treated as though the 
blades extended to the shroud because of continuity con- 
siderations. 





Distance along normal to streamline, n, FF 


(a) Parameter a. 


- I a dn 
(b) Exponential be ~o 


(c) Welght-flow factor. 
Firaurg 3.—Plots used for numerical integration of equations (10) and (11). 
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Fic-t RE 4.—Plots used for computing new streamlines. 


RESULTS AND DISCUSSION 


A solution of the flow conditions throughout the passage 
from hub to shroud is shown in figure 5. 


Streamlines.—The flow streamlines in the meridional 
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(a) Streamlines for compressible flow. 
FIGCEE 5.—Fhow analysis for centrifugal impeller in meridional plane. 


plane are shown in figure 5 (a). Although the solution has 
been made for equal quantities of fow between streamlines, 
little information can be derived from the streamline con- 
figuration, as shown, other than the direction of flow because 
the variation in distance between streamlines is due to a` 
variation in blade thickness and radius as well as in velocity. 
Near the outlet where blade thickness and radius have little 
effect, however, the streamline spacing indicates considerable 
variation in velocity across the passage. 


Velocity ratio.— Lines of constant velocity ratio Q relative 
to the impeller are shown in figure 5 (b). Along the shroud 
(fig. 6), the velocity is near sonic at the inlet (Mach number, 
0.973) with immediate deceleration just inside the impeller 
and continued deceleration to the outlet, a condition that is 
highly conducive to boundary-layer build-up and flow separa- 
tion. A similar condition exists along the straight part of 
the hub (A to B, fig. 5 (b)). Along the curved part of the 
hub (B to C), the velocity variation is small with a slight 
acceleration existing near the impellef outlet. | 

The relative-velocity distribution as obtained in the solu- 
tion represents the mean velocity from blade to blade. 
According to the analysis of reference 1, the velocity along 
the suction surface is much larger than that along the pres- 
sure surface; the velocity varies somewhat uniformly across 
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FIGURE 6.—Concluded. Flow analysis for centrifugal impeller in meridional plane. 


the passage. If the blade unloads at the outlet, as it must 
to satisfy the Kutta condition of tangency at the tip, the 
decelerations along the suction surface of this blade are of 
necessity much greater than that shown for the average 
velocity (fig. 5(b)). Such adverse velocity gradients could 
easily cause considerable separation along the flow path. 
At the impeller outlet, these flow separations would result 
in poor flow conditions with consequent losses at the diffuser 
inlet. 


Pressure ratio.—Lines of constant pressure ratio are 
shown on figure 5 (c). The pressure lines approximately 
parallel the normal lines and indicate that static-pressure 
readings on the stationary shroud of this impeller give repre- 
sentative values for the static pressure across the passage 
from hub to shroud. In figure 7, experimental values of 
pressure ratio along the shroud are compared with the theo- 
retical values obtained in the numerical solution. Pressure 
ratios are shown for the theoretical weight flow and bracket- 
Ing experimental weight flows. The experimental curve for 
the higher weight flow represents a point very near the maxi- 
mum. flow or flow cut-off line and shows a large pressure loss 
near the inlet. The experimental curve for the lower weight 
flow represents an operating point near maximum efficiency 
and consequently approaches the theoretical conditions. 
No pressure drop occurs at the inlet, which is in agreement 


with the theoretical curve; however, the theoretical pressure 
increases more rapidly along the shroud because of ineffi- 
ciency dur, to friction and flow separation, and because of a 
velocity increase resulting from flow-area reduction due to 
boundary layer and from a slightly higher theoretical tip 
speed. 

Maximum flow.—The maximum theoretical weight flow 
across any normal can be determined by making a series of 
estimates of velocity at the hub and obtaining the corre- 
sponding weight flows. For the numerical example, the 
greatest restriction of flow occurred just inside the inlet. 
The maximum flow over & range of tip speeds was deter- 
mined by the consideration of this flow restriction. Inas- 
much as the radius of curvature was infinite for the stream- 
lines near the inlet in the complete numerical solution, it 
was taken as infinite in computing maximum weight flow 
at all speeds. A comparison of the theoretical and experi- 
mental maximum flow over the speed range is shown in 
figure 8 (a). The slightly higher theoretical values are in 
good agreement with the experimental values over a rango 
of tip speeds up to 1400 feet per second (deviation of the 
order of 3 percent). At speeds above 1400 feet per second, 
the deviation increased to a maximum value of 7 percent at 
a tip speed of 1600 feet per second, 
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Figure 6.—Comperison of relative velocity and blade angle along shroud. 
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FIGURE 7.—Comparison of theoretical pressure ratia along shroud with experimental values 
for varying weight flow. 


The high relative velocities within the impeller, which 


caused choking, resulted from the large blade angle 8 near 
the impeller inlet. A large inlet blade angle was required 
for shockless inlet conditions resulting from the considera- 
tion of the impeller mass flow and the rotational speed. 
These high relative velocities, however, can be reduced by 
an appropriate selection of blade shape to give a more rapid 
decrease in f near the impeller inlet. A comparison of £$ 
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(a) Comparison of marimum fiow capacity. 
(b) Comparison cf blade-angle decrease from Inlet to outlet at shroud. 


FIGURE §.—-Comparison of change in blade angle and maximum flow capacity of two impellers 
having sams heb-shroud profiles and blade inlet angles. 


for the example impeller A and a second impeller B (refer- 
ence 3) is shown in figure 8(b); 8 for each impeller is plotted 
against the distance ratio L along the shroud. Impeller B 
has the same hub-to-shroud profile and inlet blade angle 6 
as that of the impeller A; however, as indicated in figure 8 (b), 
the rate of decrease in blade angle for impeller B is 
much larger near the impeller inlet. The resulting higher 
flow capacity of impeller B (fig. 8 (a}) may be attributed to 
the larger rate of decrease of blade angle and consequently 
lower relative velocities in the impeller passage near the 
inlet. 
SUMMARY OF RESULTS 


A method has been developed for a compressible flow 
analysis between the hub and the shroud of mixed-flow 
impellers of arbitrary design. Axial symmetry wes assumed, 
but the forces in the meridional plane that are derived from 
tangential pressure gradients were taken into account. 
Application of the method to @ mixed-flow impeller for 
which experimental data are available yielded the following 
results: 

i. A large variation in velocity existed across the passage 
from. hub to shroud. 


2. Adverse velocity gradients occurred along both hub — 


and shroud surfaces. 
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3. The theoretical pressure ratio increased more rapidly 
than the experimental along the shroud from inlet to outlet. 
Inefficiency of the actual compression process due to flow 
separation accounts for the greater part of this deviation. 

4, The maximum flow over a range of speeds as theoreti- 
cally determined closely agreed with the experimental. 

5. For the example calculation, values of static pressure 
at any point along the shroud were representative of those 
across the flow path from hub to shroud. 


CONCLUSION 


Analysis of flow in the meridional plane of a mixed-flow 
impeller showed that the rotational forces, the blade curva- 
ture, and the hub-shroud profile can introduce adverse 
velocity gradients along the hub and shroud surfaces. 


Lewis FLIGAT PROPULSION LABORATORY 
NATIONAL ADVISORY COMMITTEE FoR AERONAUTICS 
CLEVELAND, Onto, April 13, 1960 


APPENDIX A 
SYMBOLS 
The following symbols are used in the calculations and | U actual impeller tip speed, ft/sec 
the figures: u tangential velocity relative to impeller (positive in 
i cos? 8 sin’@dr_ sin?8 . direction of rotation, fig. 2(c)), ft/sec 
g “paramere equal to tht r anor snatana | y through-flow component of velocity (fig. 2(c)), ft/sec 


b parameter equal to 2w sin 8 on sin 8 tan æ 


(20 sin atge) 

€ speed of sound for stagnation conditions, ft/sec 

F force per unit mass 

f factor accounting for area taken up by blade cross 
section 

g acceleration due to gravity, ft/sec? 

H absolute total enthalpy, Btu/lb 

h static enthalpy, Btu/Ib 

J mechanical equivalent of heat 

L ratio of distance m at any point to m at impeller ani 

m distance along streamline from inlet, ft 

N number of blades 

n distance along normal to streamline in n meridional ane 
ft 

P total pressure, lb/sq ft 

p static pressure, lb/sq ft 

Q — velocity ratio, g/c: 

q velocity relative to impeller (fig. 1(b)), ft/sec 

r radius of particle from axis of impeller rotation, ft 

re radius of curvature of streamlines in meridional plane 
(positive in fig. 2(a)), it | 

r’ slope of streamline in meridional plane with respect to 
axis of rotation : 

r’’ change in slope 7’ with respect to axis of oaot 

t time, sec 


ODA @ N 


JF total compressor flow rate, lb/sec 
flow rate, Ib/sec 
distance in axial direction, ft 
angle between tangent to streamline in meridional 
plane and axis of rotation, deg 
8 blade angle (fig. 1, negative for backward-swept blades), 
deg 
ratio of specific heats 
ratio of total pressure at inlet to standard sea-level 
pressure 
8 ratio of inlet stagnation temperature to standard Sea- 
level temperature 
A prerotation term equal to re(wrs+ u) 
p mass density, lb-sec?/{t* 
pe total mass density, lb-sec?/ft* 
ge blade angle in axial-tangential direction (negative for 
backward-swept blades, fig. 2(¢)), deg 
w angular velocity of impeller, radians/sec 


Rue 


o& 2 


Subscripts: 

a average — 

d pressure or driving face of blade 

i inlet conditions 

n indicates position along normal from hub to shroud or 


normal component 
radial component 
suction or trailing face of blade 
axial component 
tangential component 
,1 points at succeeding streamlines along normal 


APPENDIX B 
DEVELOPMENT OF FORCE EQUATION IN HUB-TO-SHROUD PLANE FOR FLOW WHEN AXIAL SYMMETRY IS ASSUMED 


Lorenz’s equations of motion relative to the rotating 
impeller, in which axial symmetry is assumed (reference 3), 


are 
1p d'r (wr+u) 





Fy or dt r 
_1dp dz 

= p 52" dB 

p du , udr 


=u rual. oF 


where the blade force F is the force per unit mass. 
From figure 2(a), 


—_* F, =F, sin a— F, cos a 
where Fa may be evaluated from Fy (appendix C) as 
_1op . d?z lò dêr 
F sm at sin a— oe COS a— yr; COS a+- 


(B1) 
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dr : . | 











dr dæ, . dv i 
Jp e Cos raa TE (Bla) 
dz . dü 
i smn æ di EOE ae. (Bib) 
dp Opdz, opdr Op _ Op | 
i dzdn Orda” Oe er vie) 
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Substituting equations (Bla) to (B1c) in equation (BI) gives 





, id d ru)? 
F are ee oat core). : ) COS a (B2) 
From figure 2(a), 
Jant dt 
Q == F. 

Equation (B2) therefore reduces to 

3 dp ae (wr+u)* dr 7 (B3) 

o z + dn z 


n r: r 


APPENDIX C 


PROCEDURE FOR EVALUATING BLADE FORCE F, FOR IMPELLER OF ARBITRARY BLADE SHAPE HAVING RADIAL BLADE 
ELEMENTS 


For an impeller with radial blade elements, a true picture 
of the blade curvature lies on a circular cylinder having the 
same center line as the impeller. The blade force F normal 
to the blade (fig. 2(b}) has the component Fs, which may be 
determined from time rate of change of angular momentum 
in the direction of rotation. The component Fe as given in 
appendix B makes the angle » with the force F. 

From figure 2(¢), 








Fam — Ee tan eosin g 
and 
H tan 
tan o= = p 
p COSa COS a 
giving 


F= — F; tan £ tan a 


From appendix B, 


pyu Y dr dr 


aod ae 


Substituting dm=v dt and dr/dm=sin æ gives 
F= (t sin a+ 2w sin a) 


du u 
re (T+ 
x v tan 6 tan e 7 „S at -2u sin æ) 


: du ur . 
f,=—v tan f tan ( 2esin ata —— sin a tan a tan 6 


From figure 1, 


F,=—gqsin 8 tana (20 sin atsejt sin? 8 sin a tan a 


APPENDIX D 


METHOD OF EVALUATING BLADE SURFACE VELOCITIES 


From appendix C, 


where F; is the blade force per unit mass in the tangential 
direction. Substituting r=g cos $, u=q sin $, and dm= 


ag 


dz 
E a z giv es 
COS 


q sin 8 sin ga 
r 


F'e=q cos f (Se cos a-t +2 sin a} (D1) 
The quantity Fep is the blade force per unit volume, which 
when multiplied by the circumferential distance between the 
blades gives the blade force per unit area or the pressure 


difference between the blades; thus, 


Feo F 27r 


We (D2) 


=<Ap 


Substituting equation (B1) in equation (B2) gives 


IERES S20 sin a) 
(D3) 


apm on cos o( 2472) (BE os 


Assuming a linear variation of pressure from blade to blade 
yields 


P:S Pa a 
å 
Pe= Pat 


Dividing by the total pressure P; upstream of the impeller 
inlet gives 








Ps _ Pa AP 

PP eP. (D4) 
and 

p Pa, AP x 


Equation (11) is 


—ifar\? qV x 
p= pus $75 (2) -(£) + 
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For zero prerotation at the inlet (A=0) and with | and 
I 


the expression for velocity becomes 
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